Energy current magnification in coupled oscillator loops 



Rahul Marathe,^'B Abhishek Dhar,^'[l] and A. M. Jayannavar^'0 

' Departamento de Fisica Atomica, Molecular y Nuclear and GISC, 
Universidad Complutense de Madrid, 28040 -Madrid, Spain 
^Raman Research Institute, Bangalore 560080, India 
^Institute of Physics, Bhuhaneshwar 751005, India 
(Dated: September 23, 2010) 

Motivated by studies on current magnification in quantum mesoscopic systems we consider sound 
and heat transmission in classical models of oscillator chains. A loop of coupled oscillators is 
connected to two leads through which one can either transmit monochromatic waves or white noise 
signal from heat baths. We look for the possibility of current magnification in this system due to 
some asymmetry introduced between the two arms in the loop. We find that current magnification 
is indeed obtained for particular frequency ranges. However the integrated current shows the effect 
only in the presence of a pinning potential for the atoms in the leads. We also study the effect of 
anharmonicity on current magnification. 
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Mesoscopic systems of micron size have recently been 
studied extensively. In these systems, at low tempera- 
tures the mean free path of an electron can exceed the 
sample dimensions, thus maintaining the coherence of the 
single particle wave function across the entire sample. In 
such coherent systems several novel effects have been ob- 
served which do not have any classical analogue [l], 
An interesting example is that of a mesoscopic conduct- 
ing loop connected through metallic leads to two separate 
electron reservoirs at different chemical potentials hl and 
/Xij respectively. A current is established between these 
two reservoirs. In the presence of this transport current 
it has been shown that current magnification (CM) can 
occur in the loop U 0] ■ In this case under appropriate 
conditions depending on the Fermi energy, the current in 
one arm of the loop exceeds the total current across the 
system. For current conservation at the junctions, the 
current in the other arm flows in the reverse direction, 
i.e opposite to the transport current in the leads. The 
predicted circulating current density can be very large 
y and has been termed as giant persistent current [5|- 
Several studies have shown the existence of CM in a va- 
riety of models including multi-channel systems Q , and 
systems with spin f?], [§1 and thermoelectric Q currents. 

Motivated by these studies on quantum systems, in 
this paper we address the question as to whether CM 
can occur in classical models of energy transport in oscil- 
lator networks. Energy transport in such systems is by 
lattice vibrations which in the long wavelength limit cor- 
responds to sound-waves. Since CM is basically a wave- 
phenomena, one expects that it should be observable 
also in oscillator networks. Heat conduction in such net- 
works have earlier been studied by Eckmann and Zabey 
[lot who also discussed the possibility of circulating cur- 
rents. The system we consider is schematically shown 
in Fig. ((!]). It consists of a loop formed by particles, 
all with masses m and connected by springs with stiff- 
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FIG. 1: (color online). Models studied: (a) Oscillator loop 
connected to infinite leads formed by coupled oscillators, (b) 
System with finite leads connected to heat reservoirs. 



ness k . Two Leads [(1) and (4) in the figure] formed by 
semi-infinite oscillator chains are connected to the loop 
in an asymmetric way such that the loop has two arms 
(2) and (3) of unequal lengths. We consider two situa- 
tions. First we study the transmission of single frequency 
plane waves across this geometry. Next we consider the 
case where a band of frequencies are fed into the leads by 
connecting them to heat baths kept at different tempera- 
tures. We also look at the effect of anharmonicity on CM. 
Since CM is a wave- phenomena one expects it to be sup- 
pressed in the presence of any thermalizing mechanism 
such as arising from inelastic phonon-phonon scattering 
due to anharmoncity. In fact for macroscopic systems 
these kind of processes would lead to an Ohmic (diffu- 
sive) behavior of the wire in which case one would not 
get CM. In the electron case one important source of 
thermalization is electron-electron interactions. However 
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it is quite difficult to study its effect on CM without using 
approximation metfiods. On the other hand for the clas- 
sical oscillator case, it is easy to introduce anharmonicity 
and numerically see its effect on CM. 

We first discuss transmission of monochromatic plane 
waves across the system. The number of particles on the 
loop are L = L2 + £3 and particles in different regions 
are numbered as in Fig. ([1^). For a wave incident from 
the left side the particle displacements are given by: 

xf'^ = Re[ (e*9'i + re-*«'i) e"*"^ ] , - 00 < < 
xf^ = Re[ {h e*«'= + n e-*«'=) e'"^^ ] , < h < L2 



xf^ = Re[ ih e'"'^ + r2 e"*''^) e"*"^" ] , < < ^3 
xf^ = Re[ t e'«'^ e"'""] , < ^4 < 00, (1) 



with q £ (0, tt) and uj = 2y ^ sin(g/2), the usual disper- 
sion relation for the harmonic chain. The six unknown 
transmission amplitudes can be found out by matching 
the solutions Eqns. ([T]) at the junctions and considering 
the equations of motion of the particles at the junctions. 
Then one gets: 
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These linear equations can be explicitly solved to get all 
the unknown amplitudes. We do not give explicit ex- 
pressions here since they are quite lengthy. The instan- 
taneous energy current from site Z to / -I- 1 is the product 
of the velocity of the (/ -t- 1)*'^ particle and the force on 
it due to the l^^ particle. Thus the energy current av- 
eraged over a time period between two neighboring sites 
on any of the four regions (s = 1,2,3,4) is given by: 



One then 



/(^) = -k (c^/2^) dT 
finds = {k uj sin q / 2)T'''^\uj) where the transmission 

r2|2 and T^^) = j^gp - jrap . In 



coefficients are given by T'^-* = T*-^' — 1 



r(2) = |t2p - 



Fig. ([2]) wc plot these transmission coefficients as a func- 
tion of frequency uj for a particular choice of parameters 
k, m, L2, L3 with L2 < L3. The most interesting features 
that we see are that, for certain values of the frequency, 
the transmission T^^-* [T^^^] on one of the arms of the loop 
can be negative and when this happens there is CM on 
the other arm, i.e \T^^'> /T^^^ > 1 [\T^^'> /T^^'^\ > 1]. In 
contrast, the inset of Fig. ([2]) gives results for L2 — L3 in 
which case all transmission coefficients are positive with 
magnitudes less than one. Note that current conserva- 
tion implies /(^^ = J^^) + = j(4) and hence negative 
current flow on one arm necessarily implies CM on the 
other. 

Next we ask the question as to what happens when the 
two leads are connected to white-noise heat baths at dif- 
ferent temperatures and Tr. Do we still get negative 
heat current flow in one of the arms and CM ? For sys- 
tems with diffusive heat flow this is clearly not possible. 




FIG. 2: (color online). Plot of transmission- versus- frequency 
in different regions of the network for the case of infinite leads. 
Parameters used are k — 1, m = 1, L2 = 4: and L3 — 6. Inset 
shows the same plots for the case where there is no asymmetry 
in the loop arm lengths ( L2 = I/3 = 6). 



In the nonequilibrium steady state one would expect the 
system to be in local thermal equilibrium which means 
that the local temperatures at the junctions will define 
the direction of heat flow and this will then be unique 
(high-to-low temperature). However, for ballistic heat 
flow, as in a harmonic system, a local temperature does 
not have a thermodynamic significance and there is a 
possibility that we can get CM. We again consider the 
harmonic network shown in Fig. ([TJd) with each lead con- 
sisting of a finite number of particles, say N. The full 



network of the loop and leads thus has M = 2N + L 
particles and is described by the harmonic Hamiltonian 
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X'^MX + I X^^X 



(3) 



where X — {xi,X2, M is the diagonal mass ma- 

trix and $ is the force matrix. We label the particles as 
shown in Fig. ([T}d). To model the heat baths, the par- 
ticles at the free ends of the leads have extra terms in 
their equations of motion corresponding to a Langevin 
dynamics. Thus the particle at the end of the left reser- 
voir has an extra part —jlxi +'7l in its equations of mo- 
tion while the particle at the right end has an extra part 
—IrXm + "HR- The noise terms are Gaussian with zero 
mean and variances given by the fluctuation-dissipation 
relations {ilLMt)VL,R{t')) = '2kBTL,R'yL,RSLMSit - t'). 
For the same parameter values as in Fig. ^ we com- 
puted the average heat current in the four regions. The 
steady state heat current on a bond between sites a and 



h is given by J^'*^ 



X 



Xa ]) where (...) now 



denotes a noise average. Using the methods described 
in Ref the heat current in any part of the system 
can be expressed as an integral over all frequencies of the 
transmission function, which can be written in terms of 
the Green's function matrix G{uj) = [—oj'^'M.+^ — iujT]^^ 
where F is a dissipation matrix whose only non- vanishing 
elements are the two diagonal terms jl and jr occurring 
at positions corresponding to the two sites connected to 
reservoirs. Thus, assuming that a unique nonequilibrium 
steady-state is achieved, the heat current on any part can 



be written as [11] : 



(4) 



where the transmission coefficient T^^\uj) = 
2kiuj [ Ga,iiuj)Gbsi-uj) - Ga,M{i-o)Gb,A[{-io) ] with 
a, 6 being two adjacent sites on the region of interest. 
The site h is chosen to be on the right of a so that the 
convention used is that current is from left-to-right. 

In all the calculations we set with — 10, — 1 
and N ~ 20. The other parameters of the network are 
the same as for the data in Fig. In Fig. ([3]) we plot 
the transmission coefficients in the leads and in the two 
arms. We see that for particular values of frequencies 
the value of T'-^-' (w) is negative and whenever this occurs 
we have jT'^'^V^''"'^'' I > 1 implying CM. Next we calculate 
the total heat current using Eqs. Q. From our numerical 
calculations we find /t^' = 0.7661, J^^) = 0.37968, J^^) = 
0.3864. We have also verified these results from direct 
nonequilibrium simulations. Thus we do not get any CM. 
For several other choices of parameter values we find the 
same situation. Introducing impurities in the loop can 
enhance the asymmetry in the system. With this we find 
that while one gets CM in specific frequency ranges, the 
integrated current again does not show any CM. 
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FIG. 3: (color online). Plot of transmission- versus-frequency 
in different regions of the network for the same system as in 
Fig. ((2| but with finite leads (A'' = 20) connected to Langevin 
heat reservoirs at temperatures Tl = 10 and Tr = 1. The in- 
set shows T(cj) for the case where the parameters of the leads 
are chosen such that they have a narrow frequency transmis- 
sion band (see text for the parameter values). 



Let us now consider the effect of introducing a band- 
pass filter between the loop and reservoirs such that only 
frequencies over the range where CM occurs, are allowed 
to pass. It seems plausible that introduction of such a fil- 
ter will allow us to observe CM in the integrated current. 
We now test this numerically. For this we introduce a 
harmonic pinning potential on all sites of the leads. Thus 
for particles in the leads, in addition to the interparticle 
potential k{xi — x;+i)^/2, there is also an additional on- 
site potential V{xi) = kgxf /2. For the choice of pinning 
strength k^, = 0.35 and k — 0.1 (only on the leads), only 
frequencies in the range ut w (0.6 — 0.87) can pass through 
the leads. This is the range where, as can be seen from 
the plot of T(a;) data for the unpinned system in Fig. ([3]), 
we expect maximum CM. In the inset of Fig. ([3]) we show 
plots of 7~^i), T"^^-*, T^"^-* for this system. As expected, we 
find significant transmission in a small frequency band. 
The values of the integrated currents in this case are: 
/(I) = 0.01245, = 0.06681, = -0.05436. Thus 
in this case CM is observed even for the integrated cur- 
rent. Note that there is CM on both the arms and the 
magnification factor is more than about five times. In 
Fig. ^ we show the temperature profile in this network 
with pinning. Here we define a local temperature at any 
site through its mean kinetic energy, i.e Ti — m{xf) at 
the l^^ site. We see that the temperature profile is non- 
monotonic. As mentioned earlier this is not surprising 
since the local temperature in this integrable system does 
not have the usual thermodynamic meaning. 

We now consider the effect of introducing anharmonic- 
ity in the Hamiltonian of the loop. We consider a quartic 
onsite potential of the form V{x) = Xx'^/A at all sites of 
the loop. We fix all other parameters to have the same 
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FIG. 4: (color online). Plot of the temperature profile in 
the network for the same system as considered in the inset 
of Fig. Q without anharmonicity (A — 0) and with anhar- 
moncity (A = 0.2). The temperature across the leads and the 
two arms of the loop are shown. 



values as the pinned system studied earlier. We plot in 
Fig. ([5]) the currents I*^^-* , /^^\ 7*^^^ in different regions as 
a function of the strength of anharmonicity A. We see 
that CM decreases with anharmonicity. Interestingly we 
find that for some value of A the current on the lower 
arm (I^^^) becomes exactly zero. Note also that the cur- 
rent on the upper arm seems to vanish at some value 
of A. The inset of Fig. ([5]) shows the effect of an inter- 
particle anharmonic potential of the form \{xi — a;/+i)^/4 
between particles on the loop. We see a similar reduction 
of CM as in the onsite case though there are some qual- 
itative differences. For macroscopic systems where the 
effective mean free path due to anharmonicity becomes 
much smaller than the system size it is expected that CM 
will not be observed. In Fig. (|4]) we show temperature 
profiles for a strongly anharmonic case. In this case, in 
contrast to the harmonic case, the temperature profile is 
not non-monotonic and correspondingly there is no CM. 

Conclusions: Motivated by studies of CM in quantum 
mesoscopic systems, in this paper we have studied mod- 
els of energy transmission in classical oscillator chains. 
We have considered an oscillator loop connected to two 
external leads. The system is made asymmetric by either 
making the two arm lengths of the loop different or by in- 
troducing impurities. For single-frequency sound waves 
we find that CM is obtained over particular frequency 
ranges. For the case where the network is connected to 
heat baths which send waves at all frequencies we find 
absence of CM. This is true for various parameter sets 
that we have tried, but we do not have a proof for this 
result. We find that CM in the presence of heat baths 
can be obtained if we introduce a pinning potential in the 
leads so that only a narrow band of frequencies are al- 
lowed to pass through the loop. While we have reported 
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FIG. 5: (color online). Plots of currents J^^', J^^', as a 
function of strength of onsite anharmonicity A. The inset 
shows the corresponding plots for inter-particle anharmonicty. 

results for a small loop, we have checked that for har- 
monic systems CM is obtained for much larger sizes also. 
Finally we have looked at the effect of anharmonicity on 
CM. Our simulations show that CM is reduced but not 
completely destroyed in the presence of small anharmonic 
interactions. We also find the remarkable effect that with 
appropriate choice of parameters, the current in one arm 
can be made to exactly vanish. We expect that CM in 
phononic heat transport will be observable experimen- 
tally in mesoscopic systems such as insulating nanotubes 
where the effective mean free path for inelastic phonon 
scattering is large compared to system size. 
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